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Abstract. Among primordial magnetogenesis models, inflation is a prime candidate to ex- 
plain the current existence of cosmological magnetic fields. Assuming conformal invariance to 



2 ■ be restored after inflation, their energy density decreases as radiation during the decelerating 

c/j . eras of the universe, and in particular during reheating. Without making any assumptions 

^ ■ on inflation, on the magnetogenesis mechanism and on how the reheating proceeded, we 

show that requiring large scale magnetic fields to remain subdominant after inflation gives 
^ . non-trivial constraints on both the reheating equation of state parameter and the reheating 

I energy scale. In terms of the so-called reheating parameter, we find that Ini^rad > —10.1 for 

CS ■ large scale magnetic fields of the order 5 x 10^^^ Gauss today. This bound is then compared 

to those already derived from Cosmic Microwave Background (CMB) data by assuming a 
_ specific inflationary model. Avoiding magnetic field backreaction is always complementary to 

■ CMB and can give more stringent limits on reheating for all high energy models of inflation. 

CN ■ For instance, a large fleld matter dominated reheating cannot take place at an energy scale 

lower than typically 500 GeV if the magnetic fleld strength today is -Bq = 5 x 10~^^ G, this 
scale going up to 10^° GeV if = 10"^ G. 
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1 Introduction 

The origin of large scale magnetic fields in the universe is still a mystery [1]. There are 
evidences of their presence in the intergalactic medium and this strongly suggests that the 
origin of these fields might be primordial [2, 3]. Magnetic fields in galaxies have a strength 
ranging from 5 to 100 /iG [4] while the strength detected within clusters of galaxies is of the 
order 1 to 10 /iG [5] . Moreover limits on magnetic fields in the intergalactic medium have 
been recently derived using combined data from the Atmospheric Cherenkov Telescopes and 
the Fermi Gamma-Ray Space Telescope on the spectra of distant blazars. Assuming blazars 
are injecting both gamma and cosmic rays, Ref. [6] reports the two-sigma confidence interval 
from 1 X lO"^'^ G to 3 X 10"^^ G. Without assuming cosmic rays production, solely the lower 
limit remains model independent [6]. Other data from HESS and Fermi imposes a lower 
bound of 5 X 10~^^ G [7-9] while Faraday rotations give an upper limit as 10^^ G [10]. 

Inflation is a candidate of choice to generate primordial magnetic fields [11], because, as 
for large scale structures and Cosmic Microwave Background (CMB) anisotropics, it provides 
the means of producing effects on scales larger than Hubble radius today, starting from 
processes on much smaller length scales. As for curvature perturbations, the idea is the 
amplification of quantum vacuum fluctuations of the electromagnetic field during inflation. 
However this is not straightforward since the standard electromagnetic Lagrangian C = 
— jF^'^F^y is conformally invariant. As a result, for a comoving observer in a Friedmann- 
Lemaitre universe, mag netic fields B^' = -^'"''''^ F^fs u y always decrease as \/ a? ^ where a is 
the scale factor. In order to overcome this problem one must break the conformal invariance 
of the electromagnetic field during infiation. Several mechanisms have been proposed to do 
so [11-22]. Let us however stress that a whole class of these models has been recently shown 
to spoil perturbations theory during the radiation era [23, 24]. 

Once conformal invariance is broken, and the field is amplified, another issue that one 
should consider is the backreaction of the generated field onto the background. Already 
during inflation, backreaction can be problematic as the electromagnetic stress can overcome 
the inflaton energy density and this gives very strong bounds on the maximum strength 
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of the primordial fields [25-30]. Moreover, the time evolution of the electromagnetic stress 
during inflation depends on the way conformal invariance is broken thereby rendering the 
backreaction problem model dependent. 

Nevertheless, after the end of inflation, mostly all primordial magnetogenesis models 
assume that conformal invariance is restored such that, without source, the magnetic field 
decreases inversely proportional to the scale factor squared. This will be our working hy- 
pothesis in the following. Let us however mention that magnetic fields could also be affected 
by other mechanisms occurring after inflation [31-33], but these are expected to affect length 
scales not much larger than the Hubble radius at that time. As a result, focusing on super- 
Hubble modes, their associated energy density Pb = -B^/2 still behaves as radiation and the 
strength of the magnetic field today is simply redshifted since the end of inflation 



Here Bq and -Bend E'-re respectively the magnetic field on Hubble scale today and at the end 
of inflation while Zcnd = oo/'^cnd — 1 is the redshift at which inflation ended. In this paper, 
we point out that Zcnd depends on the properties of reheating [34-38] and therefore the value 
of magnetic fields today is connected to the reheating epoch for all models of inflationary 
magnetogenesis. In particular we show that requiring large scale magnetic fields not to have 
backreaction after the end of inflation, yields a lower bound on the reheating parameter 
Ini^rad > —10.1 for i?o = 5 X 10~^^ G. The reheating parameter has been introduced in 
Refs. [39-41] and can be expressed under various equivalent forms such as 



Here "reh" and "end" stand respectively for the end of reheating and inflation while u^reh is 
the mean equation of state parameter during reheating. The quantity AA^ = N^^h — ^end = 
ln(a].ch/ocnd) is the number of e-folds reheating lasted. Requiring to be consistent with 
standard cosmology, namely that reheating occurs before Big-Bang Nucleosynthesis (BBN) 
and after inflation, one gets pnuc = (lOMeV)'^ < /?reh < Pcnd < (10~^-^p;)^- The upper 
limit comes from the observed amplitude of the CMB anisotropics^. Moreover, the positivity 
energy conditions in General Relativity impose that — 1/3 < w^ch < 1 such that the reheating 
parameter could take any value in the range Ini^j-ad ^ [~35, 12]. Our result, IniZj-ad > —10.1, 
is therefore a non-trivial lower bound. As Eq. (1.2) emphasizes, this limit can be propagated 
either into a lower bound on the energy density of reheating, or number of e-folds, if one 
assumes the equation of state to be known; or the converse. 

Let us mention that the seven year Wilkinson Microwave Anisotropics Probe (WMAP7) 
CMB data have been shown to constrain the reheating parameter, but provided an inflation- 
ary model is specified [41-47]. Our result does not need this assumption and will be compared 
to CMB in Sec. 3. 

Equation (1.1) is valid at any length scale not affected by post-inflationary mechanisms 
and more recent astrophysical processes. Consequently, by considering magnetic fields on 
Hubble length scales today, our bound is a necessary condition for avoiding backreaction. 
For definite primordial magnetogenesis models, the primordial magnetic field spectrum is 
known such that tighter constraints may be derived [20, 48, 49]. Up to our knowledge, the 

^Mp, = I/^/SttGn is the reduced Planck mass. 




(1.2) 
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link between magnetic field backreaction and reheating duration has only been discussed in 
Ref. [20] in the context of scale invariant models of inflation breaking conformal invariance. 
In this work, by assuming a magnetogenesis model, the authors use the reheating bounds 
coming from CMB to get information on the allowed values of the primordial magnetic field 
spectrum, its spectral index and the energy scale of inflation. Here, we do not assume any 
magnetogenesis model and consider the magnetic backreaction problem in an inverted way 
compared to Ref. [20]. By preventing the Hubble mode to backreact after inflation, we 
extract some information on the reheating in a model independent way. As discussed in the 
conclusion, our results can also be applied to non-inflationary magnetogenesis model provided 
some parameters are reinterpreted according to the context (see Sec. 4). 

The paper is organized as follows. In Sec. 2, we discuss the origin of the i?rad constraint 
and derive it in terms of Bq. In Sec. 3, its implications on the energy scale of reheating are 
explored by making some extra assumptions on the equation of state parameter w^ch- We 
finally compare it with the current WMAP7 constraints on i?rad for some specific inflationary 
potentials and conclude in Sec. 4. 

2 Reheating and magnetic fields backreaction 

We assume that magnetic fields are created (for instance by quantum vacuum fluctuations) 
and amplified during inflation by some mechanism that we let unspecified. Our working 
hypothesis is however that conformal invariance is restored at the end of inflation [11-19, 21] 
such that B decays subsequently as 1/a^. Therefore the strength of the field today Bq on 
large scale is given by Eq. (1.1). The actual value of oo/aend = 1 + -^end depends on what 
happens in the universe after inflation and it is therefore related to the properties of the 
reheating period. 

2.1 Reheating parameter 

Assuming instantaneous transitions between inflation, reheating, radiation and matter era, 
the reheating parameter R^ad is defined by [39] 

Rr.d^^(^Y\ (2.1) 

where "reh" means at the end of reheating, which is also the beginning of the radiation era. 
With this definition, i?rad quantifies the deviation the reheating may have compared to a 
pure radiation era (i?rad = 1 in that latter case, or if reheating is instantaneous). From this 
definition, one can immediately evaluate the redshift at which inflation ended 

l + .end = (l + W^^ = 7?^ ^ > (2.2) 

where p.y is the energy density of radiation today^. Introducing the instantaneous equation 
of state parameter during reheating, w = P/p, using only energy conservation one has [41] 

Preh = Pendexp (-3 Z"^"' [1 + w{a)] d(lna)| = Pende-''^'^('+^-'^\ (2.3) 

^We have neglected a small correction eventually coming from non-relativistic neutrinos today and — 
3-ffo Afp, f^rad where firad the density parameter of radiation today (Slrad — 2.5 x 10~^h~^). 
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where AA^ = Nj-^h — N^nd ^ find Wj-ah is the mean equation of state parameter defined as 

w^reh = / w{a)d{lna). (2.4) 

-'Vend 

Plugging Eq. (2.3) into Eq. (2.1) gives back the two equivalent forms of Eq. (1.2). 
2.2 Avoiding magnetic fields backreaction 

After the end of inflation the energy density of the produced magnetic field Pb scales as 
radiation. In order to avoid backreaction on the background, we should consider two cases. 

A possibility is that the reheating era has W^ch ^ 1/3, and the energy density of the 
universe p during reheating decays faster than radiation such that backreaction on the length 
scales of interest is avoided for 

^^^ = ^<1, (2.5) 

Prch P'y 

which is a trivial statement. In other words, magnetic fields would have an energy density 
higher than photons today. We can nevertheless convert this bound into magnetic field values. 
Using Planck units (with /io = 1), one has 1 G ~ 3.3 x 10"^'^ M^^ and p^ ~ 5.7 x lO'^^^M^^ ~ 
5.2 X 10~^^G^. Any (homogeneous) magnetic field higher than 3/LiG would then gravitate 
more than photons today, but also at any time during the radiation era. As can be checked in 
Eq. (1.2), Wreh ^1/3 implies Rmd ^ 1 and all these reheating models are thus not constrained 
by magnetic field backreaction. 

The other possibility is Wrch < 1/3 for which, moving back in time, suggests that the 
energy density of the magnetic field Pb may dominate over p during reheating. The non-trivial 
request that there is no backreaction of the generated magnetic field onto the background, 
namely that the energy density Pb{z) remains smaller than the background energy density 
p{z) at any time, can be summarized as 

/'Bend < /'end • (2.6) 

A magnetic field today of value Bq corresponds to an energy density at the end of inflation 
given by 



1 Bq Pend 

2<d Pi 



Pb.... = -.^'-^ ■ (2.7) 



From Eqs. (2.6) and (2.7), we immediately get the bound on the reheating parameter: 

(2.8) 

Let us notice that Eq. (2.8) is totally independent on the model of inflation and only requires 
conformal invariance to be satisfled during the decelerating eras. Taking the intergalactic 
measurements as flducial values for the large scales magnetic fleld today, Eq. (2.8) finally 
gives 

So = 5xlO~^^G ^ IniJrad > -10.1, 

(2 9) 

Bo = W~^G =^ IniJrad > -4.0. 
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2.3 Reheating energy scale 

We can now plug Eq. (2.8) into Eq. (1.2) to get some information on the energy scale and/or 
duration and/or the equation of state parameter of reheating. For instance, preh can be 
expressed in terms of Pendi ^^reh and Rrad as 



Preh — ^end(-Rrad) 



such that the above limit on i?rad provides 



Preh > Pond 



(2p,)V2 



''(l + ™rch) 



(2.10) 



(2.11) 



Thus, if large scale magnetic fields have an inflationary origin, this formula tells us that 
once the energy scale of inflation is known, we have a lower bound on the energy scale of 
reheating. From Eq. (2.11), it is clear that if inflation occurs at very low energy scale, the 
lower bound on prch could be weaker than what BBN already tells us (prch > Pnuc)- Notice 
that the power goes to infinity for uJreh — ^ 1/3 such that, at fixed Bq, the lower bound goes 
to zero. This is expected since, as already mentioned, there is no non-trivial constraints from 
magnetic backreaction if reheating is radiation dominated. On the other hand, at fixed Pend) 
Eq. (2.11) shows that the lower limit on prch can be very sensitive to the value of Bq, unless 



-1/3 (almost infiation). 



In the next section, we compare these results to those derived using CMB data and 
show that they can be tighter for some models and always complementary. 



3 Comparison with CMB bounds 

In infiationary cosmology, CMB anisotropics directly probe the primordial perturbations 
generated during infiation. In a given infiationary model, various CMB observables such 
as the spectral index, its running or the tensor-to-scalar ratio depend on the number of 
e-folds before the end of infiation at which observable wavenumbers crossed the Hubble 
radius. This number depends on the post-inflationary evolution of the universe, and thus on 
reheating [37, 50]. Using the WMAP seven years data, it has been shown in Ref. [41] that 
reheating is actually constrained for two classes of infiationary models, the so-called large 
field and small field models. For our purpose, it is important to stress that, as discussed at 
length in Ref. [39], CMB ends up being sensitive to the rescaled parameter 

1/4 

R = Rra.dY7~^ (3-1) 

which involves an extra factor compared to the reheating parameter. 

3.1 Large field models 
3.1.1 Generic reheating 

In large field models, the infiationary potential has the following form 
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where M is the energy scale which fixes the amphtude of the CMB anisotropies and p is 
a free index. For large field models, both R and pcnd have been shown to be constrained, 
and after marginalization over 0.2 < p < 5 and over the standard cosmological parameters, 
Ref. [41] reports the two-sigma confidence intervals from WMAP data: 

In R > -28.9, 4 x 10^^ GeV < p^^^ < 1.2 x 10^*^ GeV. (3.3) 

From this equation, taking the central value ~ 8 x 10^^ GeV we get from Eq. (3.1) that 
IniJrad > —23. This has to be compared to the magnetic field limit of Eq. (2.9). As a result, 
reheating in large field models is currently more constrained by the magnetic fields bound. 

3.1.2 Assuming an equation of state 

Although Eq. (3.3) does not make any assumptions on the reheating, large field models are 
expected to end with parametric oscillations around the minimum of the potential. In that 
case, the equation of state parameter during reheating reads uJreh = (p — 2)/(p+2) [34, 35, 37, 
51]. Let us derive the bound on the energy scale of reheating given by infiationary magnetic 
fields for different choices of uJreh- 

One can argue that, by Taylor expanding the potential, and for small enough field 
values, any potential should behave as t/)^. In that case, parametric oscillations proceed with 
Wrch — and the reheating era expands as a matter era [34]. Using CMB to get perid given 
above, assuming w^ch = 0, the magnetic bound of Eq. (2.11) gives 

So = 5 X 10"^^ G pj/^ > 490 GeV, 

Bo = 10-1^ G ^ pI^^ > 1.4 X 10^ GeV. (3.4) 
Bo = 10-^ G ^ pj/^ > 4.3 X 10^° GeV. 

Another possible choice is to assume that symmetries ensure that the potential (3.2) is 
valid at all field values, then one has ItJreh = (p — 2)/(p + 2) for all p. Taking a negative W^ch 
actually yields stronger bounds. For instance, with the quite extreme value Wrch = —0.3, we 
find 

Bo = 5x 10"^^ G ^ pI!^^>1.1x 10^^ GeV, 

5o = 10-14 G ^ p^Yh > 1.6 X 10" GeV, (3.5) 

Bo = 10-9 G ^ pI!^^ > 9.1 X 10^3 GeV. 

As already noticed, the sensitivity to the Bo values goes down when uJrch approaches accel- 
eration (Wreh — ^ -1/3)- 

The numbers of Eq. (3.5) cannot be straightforwardly compared to those coming from 
CMB alone as Ref. [41] gives only marginalized results over all values of p. In order to be 
consistent, we have reproduced the same CMB analysis as in Ref. [41] fixing the power p of the 
potential at either p = 2 or p > 1 (corresponding to w^ch = or Wrch ^ —0.3, respectively). 
For the sake of clarity, we give only the two numbers we are interested in^. For p = 2, and 
thus Wreh = 0, we find the two-sigma CMB lower bound p^^^ > 70 GeV, which ends up being 
less stringent than all magnetic bounds of Eq. (3.4). Doing the same for p ^ 1 and therefore 
w^rch ^ —0.3, CMB gives p^^^^ > 10^ GeV which is again less stringent than Eq. (3.5). 

^The CMB bounds on pend remain typically the same. 



-6- 



3.2 Small field models 



Small field models have a potential given by 



Vi(p) = 1 




(3.6) 



and the CMB data analysis with WMAP7 data gives the lower two-sigma limit [41] 



Interestingly, for a generic reheating, there is no CMB bound on pcnd for small field models 
such that i?rad alone remains unconstrained. As a consequence, Eq. (2.8) is complemen- 
tary since it gives some information on i?rad- Erom Eq. (3.1), one has Ini? = IniJ^ad + 



if pj^^ > 6 X 10^^ GeV (for fiducial i?o = 5 x 10~^^G). If the energy scale of inflation is 
lower than this value, CMB remains more restrictive than avoiding magnetic backreaction, 
as expected from the previous discussions. 

Using the highest possible value of Bq = 10"^ G, the magnetic bound becomes stronger 
than CMB limits only if p^j^ > 1.3 x 10^'^ GeV. Conversely, there is no magnetic field 
backreaction problem on Hubble scale today in all small field models fitting CMB and having 
an energy scale lower than this value. In that case, Eq. (2.8) is automatically verified. 

4 Conclusion 

We have shown in this paper that avoiding magnetic field backreaction yields some non-trivial 
constraints on the reheating epoch. Our results have been derived under the assumption 
that conformal invariance is restored after inflation and can be summarized in terms of the 
reheating parameter iJj-ad as Ini^^ad > —10.1 for Bq = 5 x 10~^^ G [see Eq. (2.8)]. 

Compared and combined with CMB data, we have shown that these conditions translate 
into a lower bound for the large field reheating energy scale when one assumes Wreh to be 
known. For any large field model having a matter dominated reheating, we find the lower 
limit for p^^^ ranging from 490 GeV to 4.3 x 10^*^ GeV for magnetic field values varying from 
5o = 5 X 10"^^ G to 10"^ G. For the small field models, we find that avoiding magnetic 
backreaction is more constraining than CMB only if the energy scale of inflation remains 
large enough, typically higher than 10^^ GeV. In fact, as suggested by Eq. (1.1), magnetic 
field backreaction is all the more important for Zg^d large. For this reason, the magnetic 
bounds are expected to be important on all high energy inflationary models while being easily 
satisfied for those at low energy. Let us stress again that, although model independent, our 
bound is a necessary condition in the sense that it must be always satisfied. According to the 
shape of the primordial magnetic field spectrum, backreaction may eventually be stronger on 
smaller length scales. In this case, if the spectrum is known, Eq. (2.6) still applies provided 
one uses the magnetic field value at the scale the spectrum is maximal. 

Finally, the above analysis can be straightforwardly generalized to any non-conventional 
post-inflationary thermal history. As discussed in Ref . [45] , if the universe evolution incorpo- 
rates a new X-era, in addition to the reheating, one can deflne a new parameter exactly as 
in Eq. (2.1). All of our previous results would still apply to the combination RradRx instead 
of ^rad- In fact, Eq. (2.8) is equally applicable to any model of primordial magnetogene- 
sis having a decelerating era preserving conformal invariance and occurring before radiation 



lnR> -23.1. 



(3.7) 
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domination. In that case, pcnd has to be understood as the energy density of the universe 
when this era starts and p^ch when it ends. 
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